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ABSTRACT. The core of quantum tomography is the possibility of writing a generally 
unbounded complex operator in form of an expansion over operators that are generally 
nonlinear functions of a generally continuous set of spectral densities — the so-called quo- 
rum of observables. The expansion is generally non unique, the non unicity allowing 
further optimization for given criteria. The mathematical problem of tomography is thus 
the classification of all such operator expansions for given (suitably closed) linear spaces 
of unbounded operators — e.g. Banach spaces of operators with an appropriate norm. Such 
problem is a difficult one, and remains still open, involving the theory of general basis in 
Banach spaces, a still unfinished chapter of analysis. In this paper we present new non- 
trivial operator expansions for the quorum of quadratures of the harmonic oscillator, and 
introduce a first very preliminary general framework to generate new expansions based on 
the Kolmogorov construction. The material presented in this paper is intended to be help- 
ful for the solution of the general problem of quantum tomography in infinite dimensions, 
which corresponds to provide a coherent mathematical framework for operator expansions 
over functions of a continuous set of spectral densities. 
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1. Introduction 

The state of a physical system is the mathematical description that provides a complete informa- 
tion on the system. Its knowledge is equivalent to know the result of any possible measurement on the 
system. In classical mechanics it is always possible, at least in principle, to devise a procedure made 
of multiple measurements which fully recovers the state of a single system. In quantum mechanics, 
on the contrary, there is no way, not even in principle, to infer the quantum state of a single system 
without having some prior knowledge on it M\. It is however possible to estimate the quantum state 
of a system when many identical copies are available prepared in the same state, so that a different 
measurement can be performed on each copy. Such a procedure is called quantum tomography. 

The problem of finding a strategy for determining the state of a system from multiple copies dates 
back to 1957 by Fano [0], who called quorum a set of observables sufficient for a complete determi- 
nation of the density matrix. However, quantum tomography entered the realm of experiments more 
recently, with the pioneering experiments by Raymer's group [Q] in the domain of quantum optics. In 
quantum optics, in fact, using a balanced homodyne detector one has the unique opportunity of mea- 
suring all possible linear combinations of position and momentum — the so-called quadratures — of 
the harmonic oscillator representing a single mode of the radiation field. 

The first technique to reconstruct the density matrix from homodyne measurements — so called 
homodyne tomography — originated from the observation by Vogel and Risken |Q| that the collec- 
tion of probability distributions achieved by homodyne detection is just the Radon transform of the 
Wigner function W. Therefore, similarly to classical imaging, one can obtain W by inverting the 
Radon transform, and then from W one can recover the matrix elements of the density operator. This 
original method, however, works well only in a semi-classical regime, whereas generally for small 
photon numbers it is affected by an unknown bias caused by the smoothing procedure needed for the 
analytical inversion of the Radon transform. The solution to the problem is to bypass the evaluation 
of the Wigner function, and to evaluate the matrix elements of the density operator by simply averag- 
ing suitable functions of homodyne data: this is the basis of the first unbiased topographic technique 
presented in Ref. Clearly the state is perfectly recovered in principle only in the limit of infin- 
itely many measurements: however, for finitely many measurements one can estimate the statistical 
error affecting each matrix element. For infinite dimensions there is the further problem that the 
propagation of statistical errors of the density matrix elements make them useless for estimating the 
ensemble average of some operators (e. g. unbounded), and a method for estimating the ensemble 
average is needed, which bypasses the evaluation of the density matrix itself, as was first suggested 
in Ref. |||. For a brief historical excursus on quantum tomography, along with a review on the 
generalization to any number of radiation modes, arbitrary quantum systems, noise deconvolution, 
adaptive methods, and maximum-likelihood strategies the reader is addressed to Ref. [pp. 

The most comprehensive theoretical approach to quantum tomography uses the concept of frame 
of observables, i. e. a set of observables spanning the linear space of operators, from which one 
derives contextually the quorum of observables and the estimation rule. The ensemble average (H) 
of any arbitrary operator Bona Hilbert space H is estimated using measurement outcomes of the 
quorum {Oi} upon expanding H over a set of functions fn(Oi) of the observables {Oi}. What 
makes the general theory nontrivial in infinite dimensions is the crucial role of the non linear func- 
tions f n (Oi) in making the infinite expansion convergent. Let's denote by Pj := f n (Oi),j = (n,l), 
such complete set of operators. Once you have the Pj , then the problem is reduced to the linear prob- 
lem of expanding an operator as H = Ysj(Qj> H)Pj, for a suitable "dual" set of operators {Qj}- 
(Notice that generally the index I is continuous, whence also the operator expansion.) The scalar 
product in the expansion is generally not simply the Frobenius's when we need to expand operators 
that are not Hilbert-Schmidt. The mathematical theory of frames [^, ^, [uj, |ll]] is the perfect tool for 
establishing completeness of {Pj} and for finding dual sets {Qj}- In most practical situations the 
set {Pj} is over-complete, and there are many alternate dual sets {Qj}, the non unicity providing 
room for optimization. A general theory should also classify the operators H having bounded scalar 
product with {Qi} and expansion H = YljiQjt H)Pj converging in average for a given class of 
quantum states. In infinite dimensions — e. g. for homodyne tomography — the easy known approach 
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works only for Hilbert-Schmidt operators, whence including trace-class ones — e. g. for estimating 
the matrix elements of the density operator over an orthonormal basis. For unbounded operators, 
however, such operator expansion becomes an infinite sum of unbounded terms. On the other hand, 
converging (and even finite) alternate expansions are known to exist for various unbounded operators 
[|l2|]. As we will see in this paper, the mechanism allowing "renormalization" of the expansion relies 
on the existence of null-estimators — namely operator-valued functions that have zero mean over the 
quorum — their existence being related to a group of symmetries of the quorum. The notion of null- 
estimator was first introduced in Ref. Jl3|], in the context of of homodyne tomography, where the 
quorum is made of the quadratures of the field mode — the quantum analogue of the Radon transform. 
Here the symmetry group of the quorum is the group (7(1) of rotations of the quadrature phase. The 
existence of null-estimators leads to infinitely many alternate expansions of the same operator over 
the quorum, allowing cancellations of the infinities in the expansion — a kind of "renormalization" 
procedure. 

The problem of classifying all operator expansions for a given quorum in infinite dimensions for 
given spaces of unbounded operators is a difficult one, and remains still open. It involves the theory 
of frames or even more general notions of basis in Banach spaces ^ [], [l(], a still unfinished 
chapter of analysis. In this paper we present new nontrivial operator expansions for the quorum 
of quadratures of the harmonic oscillator, and introduce a first preliminary general framework to 
generate and classify new expansions, based on the Kolmogorov construction. We hope that the 
material presented here will open the way to the solution of the problem of quantum tomography 
in infinite dimensions, leading to a general mathematical framework for operator expansions over 
functions of a set of spectral densities. 



2. Quantum tomography and quorum of observables 

The general idea of quantum tomography is that there is a set of observables {X^}^ e x on the 
Hilbert space of the system H — called "quorum" — by which one can estimate any desired ensemble 
average by measuring the observables of the quorum, each at the time, in a scheme of a repeated 
measurements. The observables of the quorum must be independent each other, namely they are not 
commuting: 

[X ( ,,X ( „]=0 ^ £'=£"• (1) 

Generally the set X parameterizing the quorum is infinite, and most commonly, is a continuum. In 
these cases, since clearly one can measure only a finite number of observables, these are randomly 
picked out according to a given probability measure on X, which, therefore, must be a probability 
space. In the following, for simplicity, we will also assume a probability density over X and denote 
it with the symbol d It follows that the ensemble average of a (generally complex) operator is 
written in the form of double expectation 

{x) = Jdn(£)(MXe\X)) t (2) 

where the generally nonlinear function f^(x\X) of the variable x has an analytic form which depends 
on the particular operator X. We will call the function f^(x\X) the tomographic estimator for X 
with quorum {X^}^ e x. If we want to achieve the estimation of (X) = Tr[pX] — the expectation 
being supposedly bounded on the state p — by averaging the estimator f^(x\X) over both quorum 
and measurement outcomes with a bounded variance, we need to have the function f^(x\X) square- 
summable over x and £, more precisely 

f dfi(0 f {dEs(x))\fs(x\X)\ 2 <oo, (3) 

J X J 

where X% denotes the spectrum of X^, and dE^(x) its spectral measure. In the following, for 
simplicity, we will consider the spectrum X$ = X independent on £. Clearly, the above square- 
summability will depend again on the state p and on the operator X. 
We first want to notice two main features of estimators: 
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(1) The estimator (x\X) is generally not unique, namely there can be many different estima- 
tors for the same operator X. This is equivalent to the existence of null estimators, namely 
functions n^(x) such that 

f dp(£)f dE((x)nt(x)= 0. (4) 

Accordingly, the estimators are grouped into equivalence classes, each class corresponding 
to an operator X. For such equivalence we will use the notation ~, i. e. we will write f — g 
or / — g ~ to denote that the two estimators are equivalent, namely they differ by a null 
estimator. 

(2) For fixed x and £ the estimator f^(x\X) must be a linear functional of X, namely 

f e (x\aX + bY)=af ( {x\X)+bf ( (x\Y), fe(x\X*) = f ( (x\X)\ (5) 

Example 1 (Homodyne tomography). |Q] The quorum is given by {X^}y ^) with X$ denoting the 
quadrature X^, = \{a) e 1 ^ + ae~ l< ^), a,a^ being the annihilation and creation operators of the 
harmonic oscillator with commutator [a, a'] = 1. Estimators for the dyads \n){m\ made with the 
orthonormal basis { \n)}, n = 0, . . . , oo are 

U(x\\n)(m\) = f + ™ M^- e ^ k2 - ika: (n + d\e ikx <i>\ri) 

dk\k\e^T k2 - akx k d L d n {k 2 ) , (6) 

where L^ix) denotes the generalized Laguerre polynomials. For the unbounded operators a and 
a* a one can check that the following are unbiased estimators 

U(x\a) = 2e" l 'x, 
U(x\a f a) = 2a; 2 - §. (7) 

The problem of quantum tomography is to establish the general rule for estimation, namely 
Definition 1 (Estimation rule). Given the quorum {X^J^x find the correspondence: 

f$(x\X) X, for every operator X on H, (8) 

where we possibly mean to find the whole equivalence class of estimators f^(X^\X). 

Before solving this task, first one needs to know that the set of observables {X^j^gx is actually 
a quorum. The easiest thing to do, however, is to derive both the quorum and the estimation rule 
contextually, starting from a spanning-set of observables — shortly observable spanning-set — namely 
a set of observables { JL}^££> in terms of which we can linearly expand operators as follows 

X= f &uc w {X)F u . (9) 
Jo 

Notice that the notion of operator spanning-set used here generalizes the notion of frames for Banach 
spaces to unbounded operators (see also the following), and is generally not strictly a frame accord- 
ing to the definition of Refs. [Q, In the following throughout the paper we will always assume 
probability distributions admitting densities. Generally, the set D is unbounded, and the measure 
d lu is not normalizable, whence, as such, the expansion ^ cannot be used for quantum tomography. 
However, generally this feature is related to the redundancy of the observable spanning-set, which in- 
cludes many observables F u that are just different functions of the same observable. Then, collecting 
the observables of the spanning-set into functional equivalence classes each corresponding to an 
observable of the quorum {X^}^ e x, one can relabel the observable spanning-set as F K ^ = f K (X%) 
with writing 

x= [ d/i(0 / du( K ) c, e (x)/„(j: c ) = / dn(0fdW), (io) 

JX JR ( JX 



y(n + d)\J_ 
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where the function fc(x\X) is the integral over the observables equivalent to namely 



f e (X e \X)= dv( K ) c, ( (X)f K (X e ). (11) 

Notice that in terms of the spectral decomposition of we can write 

X= [ d/*«) / dE e (x)f ( (x\X), (12) 
Jx Jx ( 

and since this expansion is linear in the spectral measure d Ec (x), the latter can be regarded itself as 
an observable spanning-set. Indeed, upon introducing the spectral density dEc(x) = Z^, x dx, and 
the density d /z(£) = m(£) d £, and renaming £ = (£, x) and 3 = {{£,, x), x £ X^, £ £ X} Eq. ( Jl2] ) 
can be rewritten in the same form of Eq. (^), namely 

X = J dCc' c (X)Z c , (13) 

where the new expansion coefficients are now given by 

4„(J0 = m(0/e(x|JO = m(0 / di/(«) c K)C (X)/ K (x). (14) 

For homodyne tomography the above quantities are explicitly given in Table |l|. 

For F u an operator frame, the coefficients of the expansion (^|) can be rewritten in form of a 
pairing (-|-) with a dual frame G^ as follows 

CuiX) = {G U \X), (15) 
in terms of which Eq. (jl4j) rewrites in the pairing form 

c'^(X)=m(0(W s , x \X), (16) 

with dual frame 



W e , x = / dv(K)Gtefi(x). (17) 

Jfi/: 

From the last equation it follows that the estimator itself can be written using the pairing 

fe(x\x) = (w e , x \x), (18) 

or, in terms of the original observable frame 

f t (x\X)= [ dv(K)(G ( , x \X)f K (x). (19) 

J Ac 
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TABLE 1 . Table of correspondence for homodyne tomography. 
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dfi(0 f Js(dE s (x))f s (x\^- l (X)). 



2.1. Unbiasing noise. It is possible to estimate the ideal ensemble average (X) by measuring the 
quorum in the presence of instrumental noise, when the noise map jV is invertible, or, more gen- 
erally, if there exists the right inverse of ,jY . In terms of observable frames this just corresponds to 
using a different dual frame. More precisely, one has: 



(20) 



Notice that Eq. (|2fj) means that 

<X)= / o>(£) / {&E&x))j,h{x\jr-\X)\ (21) 

J X J 

where (■} = Tr[p-] denotes the ideal ensemble average, whereas (-),y/ = Tr[,yl /T (p) • ] de- 
notes the experimental ensemble average, being the predual map of j¥ (Schroedinger versus 
Heisenberg picture). This also means that for left invertible map ,jV the noisy spectral measures 
,yV(d E^(x)) are still a quorum. In terms of the pairing in Eq. ( |lq ) unbiasing the noise is equivalent 
to use the new dual frame ,yV~ lff (W^ >x ). When JV is not right-invertible one can still estimate the 
ensemble average of operators in the range of the map. Moreover, in infinite dimension, when the 
noise CP map JV is compact its inverse map is unbounded, and one generally cannot unbias the noise 
without restricting the space of reconstructed operators. Otherwise, one has a Hadamard ill-posed 
problem, for which there are biased compromises, such as putting a cutoff on the vanishing singular 
values of jV . 

Example 2 (Pauli tomography in a Pauli channel). |Q| The operators [o a /y/2] make an observable 
orthonormal basis for C® 2 . We consider now the noise described by the depolarizing Pauli channel 

JS = (l-p)jr+^J7, (22) 

where J? denotes the identity map and S? the trace map 3?(X) = /Tr(X). This noise can be 
simply unbiased via noise-map inversion: 

^ _lt = t^—J - WT^r^- (23) 

l-p 2(1 -p) 

Example 3 (Homodyne tomography with r\ < 1). |Q] The set of displacements operators D(a) := 
e aa ~ a a , a 6 C make an observable Dirac-orthonormal frame for T \/2, where 

J S ={X £J, X = :/(a,a t ):, s.t. lim /(Q,Q)e s|a|2 = 0} (24) 

OL — >00 

where : : denotes normal ordering. In the presence of noise from nonunit quantum efficiency r\, the 
unbiased reconstruction is possible for operators in T s if 7] > (2s) -1 . In fact, one uses the new 
dual: 

D{a) -> Js-^(D(a)) = riD^ 1 ' 2 a)e^ a? . (25) 

Example 4 (Homodyne tomography in Gaussian noise). |Q| As for quantum efficiency, Gaussian 
noise can be unbiased for mean thermal photon number n < s — i. One has the new dual: 

D(a) -> JS-^(D(a)) = £>(a)e" |a|2 . (26) 



3. The case of homodyne tomography 

Before addressing the general problem of deriving a general tomographic rule for unbounded 
operators, in this section we will re-derive the known pattern function of homodyne tomography in 
order to illustrate the general concepts introduced in the previous section. 
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The starting point is the observable frame {D(a)} a6 c, in terms of which the decomposition (q) 
for Hilbert- Schmidt operators rewrites as follows 

X = f — Tr[DUa)X]D(a). (27) 

JC T 

By changing to polar variables a = (— i/2)ke 1 ^ , Eq. ( p7| ) becomes 

x= r d4 r +o ° dm &lkX e _ lkX ^ 

Jo T i-oo 4 

Eq. ( ^ ) can be used only for Hilbert-Schmidt operators, for which the trace under the integral in Eq. 
(E8I) exists. In terms of the quadrature spectral measure, one has 



X = 
where 



Jo 1" J-oo 4 Jo i-oc 

= e-^ ata -D(a;)Wo,o J D t (a;)e i0ota I W ,o = ~s p V2 . ( 30 ) 



2- X*' 



P denoting the Cauchy principal value. On the other hand, in the next section we will show that for 
unbounded operators we also have the expansion 



x= rd£ 1"° dt Tr{GHt,4>)X]F(t,<b), 

Jo n J-00 

F{t,<t>)=-^e-^ t ~ t2X 4> ) \ G{t,(j>) = ^-te* [ dO\i(l - e)te i4> ) {-i9te% 
V2tv dt Jo 



(3D 



which in terms of the the quadrature spectral measure reads 

PTT 1 / roc 

X = -2.1 dE^x)rr[XW^ x }, (32) 

Jo n J-00 

where now 

/°° 1 1 2 H t' 2 f 1 

dt-^=e-i (t -* 2x) -f-teir / dem-e^e^l-idte^]. (33) 
-ec dt J lK ' 

Alternatively, as shown in Sec. ^ by means of the frame of normal-ordered moments, one has 
the expansion 

00 

X = J2 ^ n a m Tr[gl m X] 

n,m—0 

A(f> r+™ [ » (n + m\~ 1 H n+m W2a 

Jo * J-oo [ „^1 V U ) 



(34) 



where 

min (n,m) (— IV 

g n , m = ., 77 dj ... \n-j)(m-j\ ■ (35) 

The above expansions in Eqs. ^ and @ are just examples of alternate expansions which are 
equivalent for the estimation of the expectation values of (even unbounded) observables, but can 
be very different as regards the statistical noise affecting such estimation. As a matter of fact, the 
problem of classifying all possible expansions has been never solved, and, hopefully, the results of 
the present paper may suggest a unifying approach to the solution of such a difficult problem. As we 
will see in the next section, the existence of many alternate expansions is due to the symmetry of the 
quorum of quadrature operators, and the resulting properties of null estimator functions. 
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4. Calculus with null functions 

We first notice that to a null estimator function ~ it corresponds a null expansion over the 
quorum, namely 

/ dfi(0 f dE e (x)n e (x)=0<^ f &ii{Qn i (X i ) = Q. (36) 
Jx Jx Jx 

Let us recall the ordering relation [|l4|| 



{fc,0 



:o' a : s = \ — — — — r— :o' J a J : r , (37) 

f^ 3K k -mi-w v 2 y 

where {fc,i} := min(fc, Z), and s = —1,0, and 1 correspond to normal, symmetrical, and anti- 
normal ordering, respectively. We will also write the symmetrical ordering as S{a^ h a l } := 
and the normal ordering as :a' k a 1 : := :a' a 
Then we have: 

Lemma 1 (Main equivalence relation). [ fl3| ] The following equivalence relation holds 

k ±i(fe+2n+2)0 



(38) 



Proof. Since X£ = £ E?=o { k l )S{(J l a k - 1 } c ^< 2 '- fc ), one has 



^ e ± i(fc +2 + 2„)0 x fe = q yfc n q (39) 

7T 

which is equivalent to (|38|), 
Stated differently: 

Lemma 2. The following equivalence relation holds 

H k (V2x)e ±z(k+2n+2) 1' ~ 0, Vfc,n>0, (40) 
where Hk (x) denote the k-th Hermite polynomial. 

Proof. From the definition of Hermite polynomials one has 



1 1 B n 



-t 2 + A/2t(a t e i *+ae _i *) 



Then, it follows that 



(41) 



U P e = 2 :X t- 



r^e ±<(fe+2+2 " ) ^ fe (V2^) = 0, Vfc,n > 0, 
Jo 



(42) 



which is equivalent to |i"c|). 
Moreover, we also have 

Lemma 3 (Equivalence of truncated Hermite polynomials). The following equivalence relations 
hold: 

e ±in* H (i) + j Kx) „ e ±™*H 2l+n (Kx), (43) 
where we introduced the truncated Hermite polynomial 

' mm — 2m)! 
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Proof. The two identities are just the complex conjugated of each other. Therefore, it is sufficient 
to prove the identity with the plus sign. The latter is a simple consequence of identity (f38|), namely 
x k e ±i(k+2n+2)4> ^ q gy us i n g the un-truncated Hermite polynomial, we have 

i+[[n/2]] \m/r,j i „MCo \2l + n-2m 

e^H 2l+n ( K x)=e^ £ (-) ( 2 +f^ K l 

^ ml (21 + n — 2m)! (45) 

^^2l + n-2m ^i[n + 2(l-m) + 2(m — l)]cf> 

from which it follows that all terms with m > I are equivalent to zero. 

Finally, one can show the Poisson identities (whose proof can be found in the Appendix) 

Lemma 4 (Poisson identities). The following identities hold 



f{x 2 )5„(4>) * - 



7T 



+ ■ 



x/(a; )^(0) 

7T 

/n particular, we have the identity 



\ _ e —2i4> ^ — e 2i< ^ 
xe < */(x 2 e 2 '*) + xe-^/(x 2 e- 2 ^) n 



-2i<£ 1 g2i^ 



M0) - -• (47) 

7T 



5. The Kolmogorov construction 

In this Section we present the Kolmogorov construction, and its relation with the fundamental 
identity of quantum tomography. 

In the following, by I2 we denote the Hilbert space of square summable sequences of complex 
numbers, and by L2 (X) we denote the Hilbert space of square summable functions over the space X. 
For example, X = R, and L2 (X) the Hilbert space of square summable functions on the real axis, or 
X = S , and Li2(X) is the Hardy space of square-summable complex functions on the circle. In the 
following we will focus attention only to the case X = R. Consider now a complete orfhonormal set 
of functions [v„ (x)] for L2 (X) . The completeness of the set corresponds to the following distribution 
identity 

^2v n {x)*v n (y) = 5(x-y), (48) 

n 

where 5 denotes the usual Dirac-delta. Consider now a (infinite-dimensional) Hilbert space H and 
denote by [w n ] an orthonormal basis for it. The following vector 

\v(x)} = ^u„(x)K), (49) 

71 

is Dirac-normalizable, in the sense that 

(v{y)\v(x)) =5{x-y). (50) 

Consider now another (infinite-dimensional) Hilbert space K ~ H. To the orthonormal basis [u„(a;)] 
for La(X) and [z n ] for K we associate a map from the observables Ox with spectrum X on H to 
operators in B(H. H <g> K) given by 

v(X) =^„(X)® \ Zn ), (51) 

n 

where, as usual, we define the operators v n (X) in terms of the spectral resolution of X as follows 

v n {X) = j AE x {x)v n {x), (52) 

where d Ex (x) denotes the spectral measure of X. 
For X, Y £ Ox formally we write 

v{Xy V (Y) = v^X'fv^Y), (53) 
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Consider the integral kernel K(x, y), x, y g X corresponding to the positive operator K 6 B(H) 

K(x,y) = (v(x)\K\v(y)). (54) 

For any two self-adjoint operators X, Y on L/2(R), the expression K(X, Y) is well defined in the 
following sense 

K(X,Y) = v(X)Hl®K)v(Y), (55) 

where K £ B(K) is given by 

K = ^2\z n ){wn\K\w m ){z m \, (56) 

namely we can also write 

K(X,Y) =Y J M^)H^\K\w m )v m {Y) = f dE x (x) f dE x {y)K(x,y). (57) 

This is also equivalent to say that for any expansion of K(x, y) in series of products of functions 
of single variable, K(X, Y) is defined as the same expansion, ordered with the functions of X on 
the left and the functions of Y on the right. For commuting X, Y, then K(X, Y) simply represents 
the same analytic expression of K (x, y), now substituting the operators in place of the variables. As 
an example, the identity operator K = I corresponds to the Dirac-delta kernel, and for commuting 
X, Y G Ox we have 

«(X) f u(y) = S(X - Y). (58) 
By replacing now H — > H® 2 , even for non commuting X and Y, one has 

(v(Xy ® J)(J ® v(Y)) = 5(X®I-I®Y). (59) 
Moreover, similarly to Eq. (|57|), one has 

K(X®I,I®Y) = {vixy ®I)(I® 2 ®K)(I®v{Y)) 
= ^v^X)^ ®v m (Y){w n \K\w m ) 



[ dEx(x)® [ AE x {y)K{x,y). (60) 
J X J X 



The fundamental identities of quantum tomography are obtained as an expansion of the swap operator 
E over the quorum, since for any state p and observable A one has Tr[p^4] = Tr[(p ® A)_E], where 
E\i>) ® \4>) = |0) s 

5.1. Homodyne tomography. From Eq. (p8[), it is clear the swap operator can be written as 

r dd> f +oc dk\k\ _ lkXik ikXik 
E= — I —e lfcA *(8)e *. (61) 



Then, the usual homodyne tomographic formula can be obtained by the Kolmogorov construction in 
writing the swap operator as follows 

E= f ^(v(X ) t ®/)(7® 2 ® J K-)7®u(X ), (62) 
Jo i" 

corresponding to the positive kernel 

K{x,x') = \{x\\Y\\x') = ^\k\e ik(x '- x) = -Z—L— ) (63) 

where Y is the quadrature conjugated to X, with [X, Y] = The kernel is clearly positive, since 
one has 

2 

dk, 



y]K(xi,Xj)€*€j = / 

JR 



4 1*1 



(64) 
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The tomographic formula consists in the following identity 



E= I ^-K[X^®I,I®X^,) (65) 

7T 



Using Eq. (|6C|), one can also write 

E = F— J^VniXttf ®u n {Xt) , (66) 

with 

MX<f) = 1 £ l v m {X <l ,)(n\\Y\\m) . (67) 

m 

The existence of null estimator functions can be taken into account by considering any operator N n ^ 
such that 



7T 



-v n {X^) ] (g) N n ^ = 



(68) 



and any estimation rule can be obtained by the swap operator 



with D n ^ = u n (X<t,) + N ni< j„ as follows 

Tr[pX] = f ^ V) Tripv^X^] Tr[D n .^X] 
Jo n „ 



(70) 



5.2. Spin tomography. For spin tomography the swap operator writes as follows 

£ = ^£±i I g p d^sin^e 4 ^-^^ = | ^(Ji-tJ®/,!® J2-n), (71) 

and we immediately see that the kernel can be written as follows 

K{r, s ) = l(j + l){r + J\2-e+-e-\s + J), (72) 

where r,s — — J, — J + 1 ... J, and {|n}} denote any orthonormal basis of the infinite dimensional 
Hilbert space H. Such a basis can be conveniently regarded as the Hardy space of functions on the 
unit circle, with (n\z) = z n , \z\ — 1, and 

dz 

\z)(z\ = 



2ixiz 



)< z |=jT|£|e*><e*|, (73) 



e_ denoting the shift operator e_ | n) = \n— l),ande+ = e_. By introducing the vectors \v(m)) = 
m + J), we can write 

S = / ^(^(Ji-n) t ®/)(/ H «2®A')(/®u(J2-n)), (74) 

where the operator if £ B(K) is given by 

K = (j + i)(l-C), (75) 
where C = §( e + + e_) is the cosine operator. 
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5.3. Alternate expansions. The general form of the swap operator is 

E = Y,{v(X„y <g> I)(I H ®2 ® K)(I <g> v{X v )). (76) 

Introducing an invertible operator L G B(K), we can write 

E = ^(u^)* ® J)(J H ®2 <g> 

= X] «»(^) t ®«m(X„)(n|if*i;- 1 |«(0)(«(0|i-K'*|ni), (77) 

^ ro,m,Z 

where ' s an Y orthonormal basis for K. Therefore, we have all the alternate expansions on 

the quorum 

z = Yl Tr [ i; ( x " ) f z i M; ( x - ) ' ( 78 ) 

where 

M,(x) =^u m (a ! )(*(I)|£X' i H = {z{C)\LK*\v(x)), 

(79) 

m 

6. Expansion of unbounded operators over the quadratures 

As already noticed, the swap operator in Eq . (|6l| ) provides the estimation rule just for trace-class 
operators. However, it is known since Richter Jl 2pthe following formula 



71 + 771 

n 



Jo t V2 n+m 



i<j>(m—n) ^-gQ^ 



Eq. ( J80[ ) was originally derived by using nontrivial identities involving trilinear products of Hermite 
polynomials. 

Here, we provide a much simpler derivation as follows. Using the definition of Hermite polyno- 
mials in Eq. (ft-lf), one has 

in m fn + m\ fn + m\ ffc n+m _ fc . 

° a = I n ) ;L [ k ) a ° ^ 

f n d (j) \ ^ fn + m\ n + m-k i<t>(2k-n-m) i<j>(m-n) 

k -° V I (81) 
n + m\ f n d<j> _ xn+m _ c i<)>{m.-n) 



71 + 771 



n I Jn IT 

1 

n + rn\ H d 



n I J n V2 n+m 
Similarly, for the symmetrical ordering, one derives the identity 

or +« mi fn + m\ In + m\ tn+m-l s 

S{a> a } = ( I > y I JS{a a }d im 

-l 



71 + 771 \ / ag e i^( m _ n ) f n + m J5{(qt e i ^)"+ ro -'( ae - i ^)'} 







!=0 



71 + 771 \ f n d(f) iM m -n) r,n + m v n + k , oox 

m j y v e 2 • (82) 
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For arbitrary ordering, using Eq. (p7[), one obtains 
{M} ,,,, 

tfe i \ km /s\3 , tfc _j 



-i 



Using Lemma 3, one has the equivalent identity 



k + l\ d(j) j^i-k) f /s\ u 1/2 



k j ^r^^wij H ^^-s^ ] ■ (84) 

In a similar way, one can derive the useful relation 



whence 

v n P # y „ sin[(<ft- y»)(n + 1)] 

= / ^0 \ ■ (00) 

J n sm((p-(p) 
Using Eq. (|82]), for the displacement operator one obtains 

oo oo n / \ 

n = ' n=Ofc=0 V / 

- / (JT+fc)i • (87) 

The last equation can be summed using the identity 



^ (n + m)\ a \ I Z J s \ 1 + 

n,m=0 v y s=0 d=0 v 7 s=0 T 

s=0 



which gives the estimation rule 



(88) 



/^(x^|d(q)) = , ■ ( 89 ) 

e ^ct + e 

Identity ( |89| ) should be compared with the equivalent estimator given in Ref . |^] . The identity in Eq. 
( p7| ) can be also derived by explicitly using the properties of null estimator functions, as shown in 
the Appendix. 

All estimation rules f^,(X^\X) given in the present section do not correspond to an expectation 
of X as in Eq. (^). However, we can suitably recover an expectation rule — which is generally not 
unique — for any observable. Consider, for example, Eq. (|87f). Using the following identity [Il5j 



[ dee n {i-e) m =B(n + i,m + i) = 7 — ' x "2"\ ni - ( 9 °) 

Jo 



i\m\ 
(n + m + 1)! ' 

along with 

d .m+n + 1 (91) 



(m + n + l)= — t 
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one obtains the integral form for the inverse binomial coefficient 



rn + n\ d 
ml dt 



1 i 



t / d0(t0)"(t(l-0)) m . (92) 
t=i Jo 



Then, the estimator (|87|) becomes 



t I d0exp(-2X^e # a*t6»)exp(2X0e^ai(l-6»)), (93) 
t=i Jo 



and for its spectral kernel one has 

a rl 

U(x\D(a)) = — x / d6exp(-2xe"*'a*8)exp(2xe- l4 'a(l -6)). (94) 
We can rewrite the estimator in form of expectation 

/ (a;|D(a)) = -£-x F df(2a;e i0 (l - 0)|e°" , V o, * o |23;e < *0) ) (95) 
d£ Jo 

where the vectors are coherent states. Eq. (^) corresponds to the functional form 

U(x\D(a)) = Tr[W^ x :D(a):], (96) 

where 

Ws x = 4-x( d6\2xe l4 'e){2xe l4 '(l-e)\. (97) 


In order to give an expectation rule corresponding to Eq. (|8l|), we can proceed as follows. From the 
integral representation Jl5j ] 

«-(*-*.) ^ ( (98) 
-oo V 71 " 

using identity (^), Eq. ( |8l| ) for normal ordering rewrites as follows 

a tm fl » = ^ e -^™+ ) *±t£dO(-iBte i +) n [-i{l-ff)te- i *} m 



7T J-oc \ / 27T 



^ r°° dt _i ( t- i2 x^) 2 



1" J-oo V27T 

x-^-t / d#Tr[a tm a n | - i6te^)(i(l - 6)te i ' l '\]e^ , 



o 



(101) 



dt 

where we have used matrix elements on coherent states 

(p\a a | a) = a p me 2VI 1 11 (100) 
Eq. ( pg| ) is equivalent to the following expansion for operators X admitting normal ordered form 

x= r d_0 r°° dtTr [ G t ( ^ )X ] F(t5(Wj 

JO T J-oo 

F(t,0) =-^ e -i (t - i2X * )2 , G(t,<t>) = £-te£ C d0\i(l -e)te i4, ){-i0te i % 
y2ix dt J 

7. Canonical dual for homodyne tomography 

The frame theory approach to quantum homodyne tomography gives further insight to the struc- 
ture of the quorum of quadrature observables. Given a set of vectors \v n ) in a Hilbert space, if 
the positive operator F = ~}2 n \ v n){v n \ is invertible, then the scalar product between two arbitrary 
vectors can be written as 

(4>\v) = ]r<V;K)Kh> , (102) 

n 

where the set of vectors {u n = F _1 |u„,}} is called "canonical dual" of the set and F is 

denoted as "frame operator". In other words, the set {|t>n}}, along with its dual {|w n }}, is a spanning 
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set for the Hilbert space, and provide a generalized resolution of the identity. In this section we show 
that the set of (generalized) projectors |a;)^>^{a;| over the quadratures give a frame when varying 
</>, and the expansion for trace-class operators in Eq. (Q) corresponds to using the canonical dual for 
the estimation rule. 

In the following, we will make extensive use of the isomorphism between the Hilbert space of the 
Hilbert-Schmidt operators A, B on H, with scalar product (A, B) — Tt[A^B], and the Hilbert space 
of bipartite vectors \A)), \B)) £ H ® H, with ((A\B)) = (A, B), and 

|A» = ^]TA„ m |n)®|7n) , (103) 

n m 

where \n) and jm) are fixed orthonormal bases for H, and A nm — (n\A\m). Notice the identities 

A® B\C)) = | ACB T )) , A® B^C)) = \ACB*)) , (104) 

wher e r and * denote transposition and complex conjugation with respect to the fixed bases in Eqs. 
©. 

By taking {n\x)o as real, in the | )) notation the spanning set \x}^,f,{x\ corresponds to the vectors 
on H a <gi Hi, of modes a and b 

\Qx)^(x\))) = e^ a " a - bH) \x) a \x) . (105) 

Notice the identities 

{{D(z)\x)q\x)q = exp(2ixlmz) 8(Rez) , (106) 

and 

e^ (afa - bi ' b) \D{z)}} = |D(2e 10 ))) . (107) 
The frame operator can be evaluated as follows 

A(j> Z" 00 , id>(u^ a—b^ b) I \ I \ ll/l -i<l>(at a-bH) 

/ dxe vy '\x)o\x)oa{x\o{x\e 

1" i-oo 

d 2 z [ d 2 w, 



F 



— I dx I ^ J ^\D(z))){{D(ze-^)\x) \x) 00 {x\ {x\D(we-^))){{D(w)\ 
r <H f°° Ax f ^1 f ^\D(ze^))}{{D(we l4, )\e 2lxlmz - 2lxlmw 5(R & z)S(Rsw) 

JO f J-oa J It J 7T 

-f-f dt\D(tte^)))((D(ite^)\ = - f ^-L\D(z)))((D(z)\ 

T JO 71 J-oo T J 7T \Z\ 

1 



7r|a — fct ' 

(108) 

where we used the eigenvalue equation (a — b^)\D(z))) = z\D(z)}}. The inverse of F is simply 
given by 

F' 1 = n\a- 6 f | = J d 2 z\z\\D(z))){{D(z)\ . (109) 
The canonical dual is obtained as follows 



= j d a z\z\\D(z)))((D(z)\x)+\x)+ 

/roo 
d 2 z\z\ |D(ze^)))e 2la:Imz <5(Rez) = / d k \k\ \D(ike l4, )))e 2ikx (HO) 
t/ — OO 

Hence, it follows that the usual Kernel operator corresponds to the canonical dual. 
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7.1. Alternate dual frames. The dual of the quadrature projectors is not unique. However, the for- 
mula of Li fllT] ] for characterize all possible alternate duals for bounded frames and discrete indexes 
cannot provide any new dual set. By denoting the frame as {|H(k, </))))} with H(:r, 4>) = 5(Xj, — x), 
such a formula can be formally written in the form 



\@(x,<t>))) = F \Z(x,m + \f(x,m- / — / dx'({E(x',<P')\F- L \S(x,<i>)))\f (* ,<{>}, 

JO 77 J-oo 

(HI) 

where {F~ 1 \cl(x, 4>)})} is the canonical dual, and {f(x, </>)} is an arbitrary Bessel set, namely 



d x 



-|/(£C, <^)| 2 < oo. 



The scalar product that appears in the integral of Eq. (Ill) writes 



d 2 z \z\ e -<ta*CRe.+a-') 5(R ^ylmC* e <»-«)(Re(*e'«-*> + 2*) S{Re{ze ^) }) 



(112) 



(113) 



/ 



dk \k\ e 2^[-cos W -«-x'] g( ksin ^ _ 

x: 

/OO 
dke 2ik[xcas(^)-x] S ^_^ 
-OO 

7T S(x — X ) 5-n (ij) — (f>) . 



(114) 



This bi-orthogonality relation implies that the formula (111) cannot reveal any new dual set. 



8. Frames of normal-ordered moments 

In this Section, by simply applying the frame theory, we recover some results of Refs. [ p^ , 
[l7| , where the set of normally ordered moments {a^a 1 } is shown to be complete, and related to a 
biorthogonal set given on the basis of Fock states. 

From the set {a' k a 1 } we immediately write the frame operator 

OO 

F= J2 \a ik a l ))({a ik a l \ . (115) 

On the Fock basis one has 



F= > — & + n) J + n)(fc + j (Z+j 

Z ✓ T7 I 1 I 



fc,Z ,n,j— 



= e ab (a t a!®fe t fe!)e at '' 1 



7'! 

atfet 



The inverse of F simply writes 



-a+6+ ( 1 1 \ -ab 



(117) 
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The dual set g k ,i can be obtained as follows 

\g k , l )) = F- 1 \<J h a l ))= J2 \n}\n l ) ( - 1 ^ r (n\(m\a t b t -^^^a^b^\a^a l l 



- E |n)W^( B + t |Hi + 5 > (m + i|i + i + ^ ( ; +i + 8)!( ' +i + 8)! 



n,m,t,j,s= 



E l n )l m ) — 1771 5n+t,fe+j+s<5 m+4 ,! + j +s — — ■ . , ■ ■ 

" — ' tlf. s\y/n\ml(n + t)\(m + 1)\ 



(-1)% 1 n V^ k (-iy (n + t-fc)! 



x - 1 1 Vl . -1 y . 1 

= E ^=f=f 1^) H -T — ^n,»+n-> 7- T , _ ■ 2^ 



- =0 V^! " ' t! (n + t-fc)! ^ j! (n + t-fc-j)! 



E 7= T =f l n )l m )^r-' ; ">fc-t'W 



n.m.t— v x y 



min (fc,Z) 



V — [ ^ = |fc-t)iz-t) 



(118) 



The dual set is then given by 

min(fc,i) t 00 . , t 

£ MF^ |t -' >< '- t| 'S^°' > -' |0)(0| °'-' 

The dual set is unique, and in fact one has the biorthogonal relation 

TtU, ; ,a f V] = S k , k ,S u ,. (120) 



9. Frame of moments versus quadrature distribution 

The frame of moments allows to recover the estimation rule for unbounded operators as an expec- 
tation rule with a dual operator of the quadrature projectors \x) < f >< i > {x\. In other words, by inspecting 
Eq. (|8l|), we would like to write an operator G{x, (f>) such that 

Tv[G\xA)a ]n a m ]= ( n + m \ H n+m (V2x) e^ m ~ n K (121) 

\ n J V2 n+m 

In this case, the operato r Gi x, cj>) is a dual of the quadrature projectors \x) t f >( p{x\, but is different 
from the canonical dual (1 1C ), which is divergent for unbounded operators. One has 



|G t ( !S ^)» = F- 1 F|G t ( a; ,0)» = F- 1 J2 loMX^VllG^fa;^))) 

k,l= 

00 



k,l=0 

(122) 



Then we obtain 



G\x,<f>) = e^ ata G t (a;,0)e-^ tlta , (123) 
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with 



oo min(fc,Z) 

E E 

k,i=o t=o 

oo min(A;,Z) 

E E 



t\VkW. 



z<i~ ife) (l\a 



/2 k + l 



-If 



% to tw(k-ty.(i-ty. 



(124) 



The dual G(x, d>) provides also new pattern functions for the matrix elements, as previously noticed 
in Refs. niA |l8[ |. For example, for the vacuum component one has 

k 



T3r[G t (x,^)|0>(0|] = £ (-ij J^ff a »(V2: 



(125) 



Notice that such pattern functions are generally no longer bounded for x — > ±oo, even for r\ > 0.5. 
One can check that the set G(x, <f>) is dual to the set of quadrature projectors as follows 



(1 x 



6(Xt-x))){{G{x,4>)\ 



- co r-\-oo rir j ± 

dx dx > ^±\\x')^{x'\)){{G{x,cj>)\5{x-x') 

30 J -ao JO 11 

OO oo 

E I at ' afc ))«ffM I = E I ^XGM I = =191. 



(126) 



10. Generating new frames 

We can gen erate different frames by changing the function of a — which gives the frame 
operator in Eq. ( 10S ). Explicitly, we have 

d 2 z 



-f(\z\)\D(z))){{D(z) 



w d(j) 



dt\t\f(\t\)\D(ite^))){{D(ite^)\ 



= j*&± j°° If j d * z j d 2 w\D(ze^))){{D(we^)\ 
xg(\z\)g*(\w\)e 2lxlmz - 2lxlmw S(Rez)5(Rew) 



d(f> 



d 2 z f d 2 w 



7T / ^ / dx I — I —\D{z)))(lD{ze- l,p )\g{\a}-b\)\x)o\x) 



TV TV 



x (x\ (x\g*(\a^ - b\)\D(we'^)))((D(w)\ 

JO ^ i-oo 



where the functions p and / are related as follows 

I t\f(\t\) = \g(t)\ 2 . 

Using the identity 

(a-tf)\A)) = (\aA))-\Aa*))) = \[a,A])), 



(127) 

(128) 
(129) 
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for reference basis such that a* = a, e. g. the photon number basis, we have 

\a - tf\ 2 \A}} = (a f - b)| [a, A])) = \[a\ [a, A]}}). (130) 
But the double commutator can be written in terms of the (dual) Lindblad super-operator 

[a\[a,A]] = -(-S?[a] +^W])A, (131) 

where L[W].A = W x AW - ^(W t WA + AW^W). Remarkably, this is exactly the dissipative 
super-operator of the displacement Gaussian noise, corresponding to a distributed loss compensated 
by a phase-insensitive amplification. 

With the aid of the following commutator rule 



T r * 
a,[a,e 



(132) 



we easily obtain 

|a f - b\ 2 \x) \x)o = \a- b*\ 2 \x)o\x) = -\d 2 x \5{x - X))) = |^[±A 2 ](X - x))), (133) 
where T denotes the Fourier transform 



Therefore, we have 



5 (|a + -b\)\x) \x) = \F[g(±\)}(X-x))), 



corresponding to the frame 



*{x,<t>) = ^J r \g{\»)]{x*-x). 

For example, if we choose g{x) = ^/f^exp (^— %-i' 2 ^ , we have 



H(x, 4>) — V27rexp 



-—{Xt-xy 

(7 



Notice that a function 



Z(x,<j>) = - x) 



(134) 

(135) 
(136) 

(137) 
(138) 



corresponds to a frame if the function h has Fourier transform which is invertible and bounded. More- 
over two functions h and h' will correspond to the same frame operator if their Fourier transform 
have the same module. More precisely, the frame operator will be given by 



- 6+1), f(t) = -\^- 1 [h](2t)\ 



(139) 



The canonical dual frame operator can be obtained by inverting the frame operator 

F-'KlxU^xl))) = J ^ j^\D(z)HD(z)\xU\xU 



= / ^ j^y|D(^)))e 2lrfm ^(Re 2 ) = f° 



dk 



\D{-ike i,l ')))e 2ikx 



f(\k\) 



(140) 



-F 



dk 



27T f(\k\/2) 



\ e - lkx *))e lkx = |/i v (x-X^))), 



where 



h v (x) = -kT 



(x). 



J(\k\/2)_ 

Hence, it follows that the usual Kernel operator corresponds to the canonical dual. 



(141) 
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11. Other frames 



Using frame calculus, it is easy to show that the following sets of operators are spanning sets 



B n ,<t, = I — 



2\ 3 1 



7T/ y/2 n n\ 



e- x lH n (V2X^) , 



(142) 



with corresponding frame operators 



7T z — ' 

n=0 

— ^ \ B n,d>)){{B n ,<t 



-|Z|" 



(143) 



where Z = a — bK 



Appendix 

11.1. Proof of Lemma 4. Consider the Poisson form of the Dirac delta for the 27r-interval 



5 2 A4>) = Urn £ e l " l e l " 



n— — oo 



In the following we will use e to mean e = 1 — Rescaling <j> by a factor 2 we obtain 

+00 +00 

7T ^ — ' 7T ^ — ' 

n— — 00 n — — 00 

which is also equivalent to the even folding relation 

*4 +) (^) = M^) + M^+t). 

since 



n — — 00 
1 +00 

A 2|n| i2n 



inch . 

e ^ + 1 



n — — 00 



On the other hand, we have the odd folding relation 

e^S^ (<P) = 6 2 „ (0) - 6 2 „ (</> + jt) 



(144) 



(145) 



(146) 



(147) 



(148) 



M0)-M0 + t) = ^ e W [e in * - (-)V**] 

n— — 00 

1 +00 

= 1 y e |2n + l| el (2n + l)^ e ^-) (0) 

7T ^ J 



n — — 00 



From the equivalence relations (|38|), we immediately derive the equivalence 



x ~ —a; 

7T 



E 



Inl 2incA J- 2fc 

e e = —a; 



1 

7T 



1 - (ee^)* +1 + 1 - (ee-^fy 



(149) 



1 - ee 2 ^ 



1 - ee - 2i ^ 



(150) 



1 2k 
— —X 
7T 



(ee 2 ^f +1 + (ee' 



ee 2i0 _ I ee -2i<t, _ I 



+ K £ (0) 



RENORMALIZED QUANTUM TOMOGRAPHY 



21 



where the distribution 



with support in <f> — gives 



l + e 2 -e(e 2 ^+e- 2 ^) 



7T 

i f dz 



K e ((j>)e 1 ' 



d<f> 



\-e A 



l ?! O 







2ty 1 + e 2 - ete^ + e~ 



(151) 



(152) 



wher e -di n) — 1 for n > and i9(n) — for n < 0, and the integral is performed on the unit circle. 
Eq. (150), which contains identity (kl) as a special case, generalizes as follows 



1 f /(ecV*) f(ex 2 e- 2 "») 



+ ■ 



+ M0)/(2 2 ) 



7T [l-e-ie- 2 ^ l-e^e 2 
where f(z) denotes any analytic function in 2. In a similar way we calculate 



(153) 



7T z — ' 



|n| i(2ra+l)4> 



1 

7T 



2j% + 1 id) 

=— a: e 



n——k—l 

1 - (ee 2!< *) fc+1 + 1 - ( ee - 2 ^) fc + 2 



1 5 
= —a; 

7T 

which generalizes as follows 

xf{x' 2 )^5- n {(f) — —x 



2fc + l i<A 

a; e v 



1 _ ee 2i0 1 I _ ee -2i<t, 



(154) 



+ 



i0 /(^i~ 2 " 2 ^"> 



1 - e -i e - 2 ^ 



+ 



e-^/(ea: 2 e- 2 ^) 
1 - e -i e 2 ^ 



+ Ke(0) 



+ ^(<j>)e i *f(x 2 ) 



(155) 



11.2. Alternative derivation of identity (K7p. By posing a = ire 1 ^, we have 



Z7T 



= /" V W [ 52 -(^-^) e!rX * + ^W-^-7r)c- irJf *] 

= r^7r [cos rX^ (0-^) + M^-0-t)) (156) 

it) 7T 

+ jsinrX^fevr^ — 0) — ^2^(0 — — 7r))] 

Now, we evaluate separately the cosine and the sine terms. In the following, we will denote ip = 
<t> — d>. The cosine term can be transformed as follows 



E 



(2k) 



2k -\s-2k 



(~) 2k 

(2k)! 



2\n\ 2inij] 



k=0 
k 



=E^ 2 "( 

n=Q 
00 

=E 



e 2m * + e 



n=0 
— 2imp 



+ e 



n= — k 
— 2inip * 



00 



(2ft) 



(157) 



2n ( 2intp 



+ e 



-2inib r- 

— O„o 



r X^, — cos(rX0) 

°° <'_ 1 i*v 2 ' s v 2fc 

/ in 2n v 2n \ "» V I ' A 4 

r ** h (2fc + 2n)! 
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On the other hand, the sine term transforms as follows 



^_^k r '2k+l j^2fc+l / 



fc=0 ^ n=-fc-l 



+ 11 i(2n+l)V 



oo / \k 2k+l Y 2k + l / fc fe + ! 

_ \ ^ ( j ^ / 2n + l i(2n+l)V , 2n-l 



-i(2n-l)V 

(2Jfe + l)! " ' 

k=0 K ' \n=0 n=l 

oo /_ \fe 2k + l y2k + l / k k 

^ e2 . + l e ^ + l)* + ^ £2 n + l e -^ +W ] (158) 

fc=o (2fc + l)! \ n=0 n=0 

oo oo /_ \fc 2fc + l y2k + l 

ST^ 2n+l ( i(2n+l)V , -i(2n+l)i/A ^ ^ 

"h [ >h ( 2fc +D ! 

oo / \fc 2fc+l v 2k+l 



. 2n+l / i(2n+l)V , -i(2n+l)iA 2n y2n I -> ^0 

^ \ +e r A ^ 2-j (2fc + 2n + l)! 



fc=0 



It is convenient now to make the following substitutions 

r = 2|a|, e _i * = i£, c i * = t^e i * e^r = 2iae 10 . (159) 

|a| |a| 

In this way, the cosine term becomes 

cos rX^irS^ (ip) 

\ ^ ( i *n2u i2nd> , / \2n —i2nd> j \ < ov i2n l — ) (2|o 

- (( eQ ) e + ( £Q ) e 5 ™° J ( 2X ^ E ( 2 fc + 



2w ^ (- ) fc (2|a|X,) 2fc (160) 
(2fc + 2n)! ' 



fc=0 



and the sine term 

i sin rX^irSl"^ (ip)te^ 



- E (-( ec 2 "+V (2 " +1) ^ + (ea) 2 " +1 e- i(2 " +1) *) (2X ) 2 " +1 £ B^T^Tr ' 

n— fc — H~ 71 j- 

(161) 

and putting the two terms together we have 

D(a)= [* ^fi{X+\D(a)), (162) 

with 

, B ^(-) fc (2|a|X ) 2fc 



/|(X^|D(a)) = E ((-ea*)"e i »* + {eaTe^ - S n0 ) (2X,) n E 

°= oo . \ OY \ 2k + n 

=EE l £ e (-« ) « + e e (-« ) -H (2fc + n)! 



fc=0 (2^ + ^)! 



n = fc = 



(163) 
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In the limit e — > 1 the last expression can be simplified using the reordering rule 

oc oc 

\ ^ \ ^ / n + fc k.n . k n-\-k ,~n r k k\ 

2__ l 2_^a ri+ 2k{z w t + z w t — d n0 z w ) 

n=0 fc = 

oc oc 

V \ v n / h„ k,h — k | k „ h,k — h c fc„ k\ 

= 2_^2—i ah + k y z w f +z w t — d hk z w ) 

fc=0 h=k 

oo h oo oo oc 

= E E ah +k {tz) h (t- 1 w) k + E E a h +k(tz) k (t- 1 w) h - J2 *2hz' 

h=0 k=0 k=0 h=k h=0 

CXD QC 



h=0 k = 

corresponding to 
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